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Project  Abstract  (from  application) 

Understanding  the  functional  connectivity  of  the  brain  is  of  substantial 
importance  for  basic  and  translational  neuroscience,  but  existing  techniques  to 
derive  this  functional  connectivity  have  important  limitations.  The  goal  of  this 
proposed  project  is  to  develop  methods  that  support  direct  and  objective 
assessment  of  functional  connectivity  in  high-resolution  electrophysiological 
brain  signals.  The  objective  is  to  provide  the  methodological  basis  for 
construction  of  the  most  detailed  functional  connectivity  maps  generated  to  date. 

Project  Summary 

Our  entry  point  to  the  project  was  a  dataset  of  electrocorticographic  (ECoG) 
activity  that  were  recorded  while  seven  people  listened  to  music,  a  complex 
auditory  stimulus.  We  set  out  to  develop  methods  that  could  accurately  establish 
those  functional  brain  connections  that  were  related  to  processing  of  the  auditory 
stimulus. 

For  a  particular  set  of  two  brain  signal  features,  functional  connectivity  can  in 
principle  be  readily  established  using  existing  techniques  such  as  Granger 
causality.  These  techniques  then  need  to  be  applied  to  all  combinations  of  brain 
signal  features.  Because  there  are  usually  thousands  of  such  combinations  of 
brain  signal  features,  this  raises  significant  problems  of  with  multiple 
comparisons.  Moreover,  extraction  of  ECoG  features  and  application  of  Granger 
causality  itself  depends  on  a  number  of  assumptions  that  are  usually  violated 
with  brain  signals.  While  all  of  these  issues  are  well  known  in  their  respective 
fields  of  statistics  and  modeling,  they  are  not  usually  considered  in  neuroscience. 
Thus,  our  initial  work  included  a  careful  analysis  of  these  assumptions,  and  the 
design  of  analytic  procedures  that  would  not  violate  these  assumptions. 

The  next  step  was  to  establish  procedures  that  can  determine  that  subset  of 
functional  connectivity  that  is  in  fact  related  to  the  task.  We  began  by  designing 
procedures  that  would  establish  such  task-related  functional  connectivity  by 
evaluating  functional  connectivity  in  relatively  short  time  windows  (e.g.,  l-sec) 
and  then  determining  the  statistical  difference  between  functional  connectivity 
metrics  during  the  task  and  during  rest,  i.e.,  a  within-subject  analysis. 


page  1/2 


Short  Term  Innovative  Research  (STIR)  Project,  ARO 


The  figure  on  the  previous  page  shows,  for  each  of  the  seven  subjects,  those 
functional  connections  in  ECoG  gamma  activity  (shown  in  colored  arrows)  that 
were  statistically  different  between  the  task  and  rest  condition.  They  also 
highlight  those  cortical  areas  that  were  separately  identified  as  being  related  to 
auditory,  face  motor,  or  hand  motor  function  (green,  red,  or  blue  shading, 
respectively).  Thus,  this  work  produced  the  first  comprehensive  analysis  of  ECoG 
functional  connectivity  related  to  auditory  processing.  While  the  results  are 
largely  consistent  with  expectations  given  by  current  understanding  in  functional 
neuroanatomy,  this  was  not  always  the  case  (e.g.,  the  highly  statistically 
significant  (blue)  connection  from  temporal  to  superior  prefrontal  cortex  in  the 
last  subject).  Careful  analyses  determined  that  the  likely  cause  for  such 
potentially  incorrect  result  is  the  fundamental  assumption  made  by  almost  all 
task- related  analyses  of  brain  signals  that,  with  the  exception  of  changes  inducted 
by  the  particular  task,  the  brain  is  in  the  same  state  during  the  task  and  during 
rest,  which  does  not  need  to  be  the  case. 

Based  on  this  consideration,  we  then  evaluated  a  completely  new  technique  to 
establish  task-related  brain  signal  activation  and  functional  connectivity.  This 
technique  tests  statistically  whether  a  particular  brain  signal  in  one  subject 
during  task  execution  is  similar  to,  or  causally  related  to,  a  particular  brain  signal 
in  a  different  subject.  Thus,  the  only  possibility  for  a  brain  signal  in  one  subject 
to  be  similar  to,  or  causally  related  to,  a  brain  signal  in  a  different  subject  is  if 
they  are  both  related  to  the  task.  This  method  had  never  before  been  applied  to 
ECoG  activity,  and  never  before  been  applied  to  functional  connectivity.  The 
attached  paper  describes  this  method  in  detail,  and  applies  it  to  determine  task- 
related  activation  of  ECoG  activity.  This  paper  is  currently  in  submission  with 
the  Journal  of  Neuroscience,  after  favorable  review  of  our  presubmission  inquiry. 

Finally,  we  applied  the  same  concept  to  Granger  causality  that  is  evaluated  across 
subjects,  i.e.,  asking  the  question  whether  one  particular  brain  signal  feature  in 
one  location  in  one  subject  is  causally  related  to  one  particular  brain  signal 
feature  in  one  location  in  a  different  subject.  To  date,  this  completely  novel 
approach  has  produced  extremely  exciting  results  that  we  are  currently  working 
on  summarizing  in  a  journal  paper.  The  figure  below  shows  these  results  that 
were  derived  from  all  seven  subjects.  Each  arrow  indicates  statistically 
significant  causal  relationship  between  two  areas  in  the  brain.  The  color/width  of 
each  arrow  is  proportional  to  the  significance.  All  of  these  functional  connections 
are  in  line  with  expectations  based  on  current  understanding  of  neurophysiology. 
However,  this  is  the  first  time  that  such  detailed  and  meaningful  results  have 
been  derived  directly  and  objectively  directly  from  brain  signals. 
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Abstract 

An  enduring  scientific  and  philosophical  question  has  been  to  what  extent  the  human  brain  functions  similarly 
across  individuals.  Whether  the  electrical  activity  of  neuronal  populations  at  a  particular  brain  location  is  similar 
across  individuals  performing  the  same  task  (e.g.,  listening  to  music),  and  how  this  activity  is  related  to  the  under¬ 
lying  anatomy,  are  largely  unknown.  Here  we  demonstrate,  using  subdural  cortical  recordings  in  humans,  that  the 
instantaneous  activity  of  both  alpha  (8-12  Hz)  and  high  gamma  (70-170  Hz)  amplitudes  was  strongly  correlated 
across  individual  subjects  at  neuroanatomical  locations  implicated  in  auditory  processing.  In  particular,  for  ECoG 
amplitudes  in  the  alpha  band,  we  identified  robust  intersubject  correlation  (ISC)  values  only  in  areas  adjacent  to 
primary  auditory  cortex,  which  has  been  known  to  receive  afferent  auditory  projections  from  the  thalamus.  In  con¬ 
trast,  for  ECoG  power  in  the  high  gamma  band,  we  identified  robust  ISC  values  not  only  in  areas  close  to  primary 
auditory  cortex,  but  also  in  auditory  association  and  premotor  cortices,  which  have  been  shown  to  be  involved 
in  higher-order  auditory  processing.  These  results  provide  the  first  neurophysiological  demonstration  that  neural 
activity  associated  with  auditory  processing  is  conserved  across  individuals  and  provide  additional  evidence  for 
the  distinct  neurophysiological  mechanisms  of  alpha  and  gamma  activity. 

Keywords:  auditory  processing,  electrocorticography  (ECoG),  thalamo-cortical  interactions,  alpha  and  high 
gamma  activity 


1.  Introduction 

Korbinian  Brodmann  published  a  map  of  the  hu¬ 
man  cortex  in  1909  (Brodmann,  1909),  which  revealed 
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4  that  the  brain’s  anatomical  organization  and  cytoarchi- 

5  tecture  share  similarities  across  humans.  Since  then, 
e  a  large  number  of  studies  have  associated  particular 
7  anatomical  areas  with  particular  functions.  For  in- 
s  stance,  neuroimaging  studies  have  demonstrated  that 

9  brain  metabolic  activity  detected  over  a  period  of  sev- 

10  eral  seconds  during  the  presentation  of  sounds  (Hejnar 

11  et  al.,  2007)  and  movies  (Hasson  et  al.,  2004)  indexes 

12  areas  that  are  similar  across  subjects.  Further,  human 

13  lesion  studies  have  shown  that  damage  to  Brodmann 
u  areas  41-42  or  Brodmann  area  17  leads  to  a  loss  of 

15  auditory  or  visual  function  (Peretz  et  al.,  1994;  Moore 

16  et  al.,  2001),  respectively.  However,  it  remains  un- 

17  known  whether  the  temporal  patterns  of  the  neural  ac- 
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tivity  associated  with  auditory  processing  are  similar  64 
across  humans,  and  how  this  activity  may  be  related  65 
to  the  underlying  anatomy.  We  addressed  these  two  ee 
questions  using  cortical  surface  recordings  in  seven  67 
human  subjects.  es 

69 

2.  Materials  and  Methods 

71 

2.7.  Subjects  and  Data  Collection  2 

73 

We  recorded  electrical  activity  from  the  surface  ^ 
of  the  brain  (electrocorticography  (ECoG))  of  seven  ^ 
epileptic  subjects  (2  men,  5  women)  who  were  listen- 
ing  to  a  complex  natural  auditory  stimulus  (the  song 
“Another  Brick  in  the  Wall  -  Part  1”  (Pink  Floyd,  7g 
Columbia  Records,  1979)).  These  subjects  underwent  7g 
temporary  implantation  of  subdural  electrode  arrays  to 
localize  eloquent  cortex  and  the  source  of  the  epileptic 
seizure  prior  to  surgical  resection.  The  implanted  elec-  g2 
trode  grids  consisted  of  platinum-iridium  electrodes  gg 
that  were  4  mm  in  diameter  (2.3  mm  exposed)  and  ^ 
spaced  with  an  inter-electrode  distance  of  0.6  or  1  cm.  gg 
The  total  numbers  of  implanted  electrodes  were  98, 

96,  83,  109,  58,  120,  and  58  for  the  different  subjects,  g7 
respectively.  Electrodes  were  implanted  on  the  left  gg 
hemisphere  for  all  subjects  (see  Fig.  1  for  electrode 
coverage).  ECoG  signals  were  digitized  at  1200  Hz, 
synchronized  with  stimulus  presentation,  and  stored  89 
using  the  BCI2000  software  platform  (Schalk  et  al.,  90 
2004).  ECoG  signals  were  recorded  while  the  subjects  gi 
were  listening  to  the  song,  which  was  3:10  min  long,  92 
digitized  at  44.1  kHz  in  waveform  audio  file  format,  93 
and  binaurally  presented  to  each  subject  using  in-ear  g4 
monitoring  earphones.  In  addition,  we  recorded  the  g5 
same  amount  of  ECoG  signals  while  subjects  were  at  gg 
rest  with  eyes  open.  We  visually  inspected  the  record-  g7 
ings  and  removed  those  electrodes  that  did  not  clearly  gg 
contain  ECoG  signals  or  had  epileptic  activity,  which  gg 
left  93,  86,  82,  104,  56,  108,  and  57  electrodes  for  the 
different  subjects.  The  ECoG  signals  from  these  elec- 

trodes  were  submitted  to  the  analyses  described  below. 

J  102 

103 

2.2.  Extraction  of  ECoG  and  Sound  Features 

To  determine  whether  the  brain  ’  s  instantaneous  neu-  ^  05 
ral  activity  at  a  particular  location  is  similar  across  in-  ioe 
dividuals,  we  first  extracted  the  time  course  of  two  107 
ECoG  processes  (i.e.,  amplitudes  in  the  alpha  (8-12  ios 
Hz)  and  high  gamma  (70-170  Hz)  frequency  bands).  109 
Activity  in  the  alpha  band  is  thought  to  represent  no 
thalamo-cortical  interactions  (Steriade  et  al.,  1990;  m 


da  Silva,  1991;  Zhang  et  al.,  2004),  whereas  activ¬ 
ity  in  the  high  gamma  band  is  believed  to  index  the 
activity  of  local  population  of  neurons  (Miller  et  al., 
2009;  Miller,  2010).  To  extract  the  time  course  of  al¬ 
pha  and  high  gamma  activity,  we  filtered  ECoG  sig¬ 
nals  from  each  electrode  at  each  specific  frequency 
band  using  an  HR  band-pass  filter,  and  removed  spatial 
noise  common  to  all  ECoG  electrodes  using  a  com¬ 
mon  average  reference  (CAR)  spatial  filter.  We  then 
computed  the  envelope  of  the  ECoG  signal  (i.e.,  the 
magnitude  of  the  analytic  signal)  in  each  frequency 
band,  low-pass  filtered  the  result  at  5  Hz  using  an  HR 
filter,  down-sampled  the  result  to  10  Hz,  and  removed 
the  linear  trends  from  the  resulting  time  series.  To  ex¬ 
tract  the  song’s  sound  intensity,  we  computed  the  av¬ 
erage  power  derived  from  non-overlapping  10  ms  seg¬ 
ments  of  the  song.  We  then  smoothed  the  sound  in¬ 
tensity  by  applying  a  low  pass  HR  filter  at  0.5  Hz  and 
down-sampling  the  result  to  10  Hz.  All  filtering  oper¬ 
ations  were  performed  forwards  and  then  backwards 
(using  Matlab’s  filtfilt  command)  to  avoid  an  in¬ 
troduction  of  a  group  delay.  The  following  analyses 
then  determined  whether  the  time  courses  of  ECoG 
alpha  and  gamma  amplitudes  were  similar  across  sub¬ 
jects. 

23.  Intersubject  Correlation  Analysis  (ISC) 

We  first  defined  the  brain  anatomy  of  each  sub¬ 
ject  using  pre-operative  magnetic  resonance  imaging 
(MRI)  scans,  and  the  location  of  the  electrodes  using 
post-operative  computer  tomography  (CT)  imaging. 
We  then  created  a  3D  surface  model  of  each  subject’s 
cortex  from  the  MRI  images,  co-registered  it  with  the 
location  of  the  electrodes  given  by  the  CT  images  us¬ 
ing  Curry  Software  (Compumedics  NeuroScan),  and 
transformed  the  3D  model  and  electrode  locations  into 
Talairach  space  (see  Fig.  1). 

We  then  related  ECoG  activity  across  subjects.  Be¬ 
cause  the  location  of  the  electrodes  usually  varies 
across  subjects,  we  identified,  for  each  electrode  loca¬ 
tion  in  each  subject,  all  electrodes  from  all  other  sub¬ 
jects  that  were  in  close  proximity  (<1  cm)  to  that  elec¬ 
trode.  (This  way,  brain  activity  from  one  individual 
was  exclusively  compared  to  brain  activity  from  other 
subjects.)  For  each  such  pair  of  electrodes,  we  com¬ 
puted  the  average  Talairach  coordinate,  which  resulted 
in  a  total  of  6852  locations  (314+117  per  subject  com¬ 
bination)  that  formed  the  basis  for  our  analyses.  For 
each  of  these  locations/electrode  pairs,  and  separately 
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for  ECoG  activity  in  the  alpha  and  high  gamma  band, 
we  computed  the  pairwise  Spearman  correlation  co¬ 
efficient  (intersubject  correlation  (ISC)  r)  and  its  sig¬ 
nificance  (i.e.,  p- value).  We  repeated  the  same  pro¬ 
cedure  for  all  locations  and  all  21  combinations  of 
subjects.  Next,  we  determined  all  locations  whose  p- 
values  were  significantly  different  than  chance  at  the 
0.001  level  (i.e.,  p  <  0.001  after  Bonferroni  correction 
for  6852  locations),  and  projected  the  negative  loga¬ 
rithm  of  the  corresponding  /7-values  (-logio(p))  on  to 
the  3D  brain  cortical  template  provided  by  the  Mon¬ 
treal  Neurological  Institute  (MNI). 

In  addition,  we  calculated  Granger  causality 
(Granger,  1969)  to  test  the  hypothesis  whether  al¬ 
pha  activity  statistically  causes  high  gamma  activity 
or  vice  versa.  In  this  analysis,  we  used  a  model  or¬ 
der  of  10,  which  was  suggested  by  the  Bayesian  in¬ 
formation  criterion  (BIC)  (Schwarz,  1978).  Finally, 
we  determined  the  lag  between  activity  in  these  two 
time  courses  by  identifying  the  time  of  the  peak  of  the 
cross-correlation  function. 
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3.  Results 


180 

181 


3.1.  Relevant  Cortical  Locations 

The  black  dots  and  colored  areas  in  Fig.  2a  give 
the  locations  and  corresponding  -logio(p)  values  ac¬ 
cumulated  across  all  subjects,  and  highlight  areas 
wherein  neural  activity  at  these  two  frequency  bands 
was  similar  across  individuals  during  the  task  (i.e.,  lis¬ 
tening  to  music,  left  panels)  as  opposed  to  the  con¬ 
trol  condition  of  rest  (i.e.,  relaxing  on  a  bed  with  eyes 
open,  right  panels).  There  were  a  total  of  277  signif¬ 
icant  locations  across  all  subjects,  the  two  frequency 
bands,  and  the  task  and  control  conditions.  Specifi¬ 
cally,  for  the  ECoG  recordings  during  the  task  con¬ 
dition,  29  locations  had  significant  alpha-ISC  values 
(mean  correlation  ±  standard  deviation  across  the  29 
locations:  r  -  0.14  ±  0.01,  max  r  =  0.18)  and  were 
primarily  located  in  areas  close  to  primary  auditory 
cortex  (Fig.  2a,  top-left).  Moreover,  246  locations  had 
significant  gamma-ISC  values  (r  =  0.19  ±  0.05,  max 
r  =  0.33)  and  were  not  only  located  in  areas  close  to 
primary  auditory  cortex  but  also  over  auditory  asso¬ 
ciation  and  premotor  areas  (Fig.  2a,  bottom-left).  In 
contrast,  for  the  ECoG  recordings  during  rest,  only  2 
(Fig.  2a,  top-right)  and  0  (Fig.  2a,  bottom-right)  loca¬ 
tions  had  significant  alpha-ISC  and  gamma-ISC  val¬ 
ues,  respectively.  The  number  of  locations  with  sig¬ 
nificant  alpha-  and  gamma-ISC  values  during  the  task 
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was  larger  than  those  during  rest  (Fisher’s  exact  prob¬ 
ability  test,  two-tail,  p- 4.51e-7  for  alpha,  /?=1.88e-75 
for  gamma).  These  results  are  summarized  in  Fig.  2b. 
While  ECoG  activity  during  the  task  was  correlated 
across  subjects  in  the  alpha  and  high  gamma  band,  this 
was  not  the  case  for  beta  (12-30  Hz,  0  locations),  low 
gamma  (40-65  Hz,  1  location),  and  very  high  gamma 
(250-350  Hz,  0  locations)  frequency  bands. 

3.2.  Temporal  Relationship 

We  then  investigated  the  temporal  relationship  be¬ 
tween  alpha  and  high  gamma  activity,  as  well  as  the 
song’s  sound  intensity.  For  this,  we  first  averaged 
the  time  courses  of  alpha  and  high  gamma  activity 
across  all  locations  where  activity  was  significantly 
correlated  across  individuals.  Then,  we  correlated  the 
three  averaged  time  courses  of  alpha  activity,  high 
gamma  activity,  and  the  song’s  sound  intensity  with 
each  other.  The  results  are  shown  in  Fig.  3  and  reveal 
a  significant  negative  correlation  between  alpha  and 
high  gamma  activity  (r  =  -0.47),  indicating  that  high 
gamma  activity  augmentation  in  auditory  cortical  ar¬ 
eas  during  auditory  processing  is  accompanied  by  al¬ 
pha  activity  suppression  (see  also  (Crone  et  al.,  2001)). 
We  also  found  significant  positive  correlation  between 
high  gamma  activity  and  sound  intensity  (r=  0.28),  in 
line  with  our  previous  findings  (Potes  et  al.,  2012),  as 
well  as  significant  negative  correlation  between  alpha 
activity  and  sound  intensity  (r  =  -0.21).  (/?<0.0001, 
Spearman’s  correlation,  n=1800).  The  scatter  plots 
shown  in  Fig.  3b  further  illustrate  these  relationships. 
Importantly,  Granger  analysis  assessing  information 
flow  revealed  that  high  gamma  activity  predicted  alpha 
activity  (p  «  le-8),  but  not  vice  versa,  and  cross¬ 
correlation  analysis  indicated  that  the  onset  of  high 
gamma  activity  preceded  alpha  activity  by  280  ms. 

3.3.  Spatial  Relationship 

To  answer  the  second  question  of  this  study,  we  in¬ 
vestigated  the  spatial  relationship  of  these  two  ECoG 
processes  (i.e.,  alpha  and  high  gamma  activity)  with 
the  underlying  auditory  anatomy.  To  do  this,  we  re¬ 
lated  those  cortical  areas  with  significant  alpha-  and 
gamma-ISC  values  to  the  cortical  and  subcortical  ar¬ 
eas  that  are  known  to  be  involved  in  auditory  process¬ 
ing.  Specifically,  previous  neurophysiological  studies 
have  shown  that  neurons  in  the  medial  geniculate  nu¬ 
cleus  of  the  thalamus  send  auditory  information  to  the 
cortex  through  projections  that  terminate  in  primary 
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auditory  cortex  (Steriade  et  al.,  1990).  In  addition,  254 
neurons  in  auditory  cortex  have  projections  back  to  the  255 
thalamus.  Theoretical  considerations  suggested  that  256 
thalamo-cortical  projections  (Fig.  4a)  contribute  to  os-  257 
cillations  at  low  frequencies  (6-12  Hz)  (Llinas  et  al.,  25s 
1999),  and  recent  experimental  evidence  in  animals  259 
is  beginning  to  confirm  the  important  role  of  thala-  260 
mic  modulation  of  cortical  activity  (Saalmann  et  al.,  261 
2012).  Our  results  identified  significant  alpha-ISC  val-  262 
ues  close  to  primary  auditory  cortex,  and  thereby  links  253 
anatomical  structures  related  to  auditory  processing  264 
to  oscillatory  electrophysiological  activity  and  its  hy-  255 
pothesized  origin  (Fig.  4b,  c).  Together  with  the  find-  266 
ings  that  high  gamma  and  alpha  activity  are  negatively  257 
correlated,  and  that  high  gamma  activity  predicts  al-  268 
pha  activity,  our  results  are  consistent  with  the  hypoth-  259 
esis  that  activity  in  the  alpha  band  reflects  interactions  270 
between  the  thalamus  and  the  cortex  (Zatorre  et  al.,  271 
2007;  Hackett,  2008;  Saalmann  et  al.,  2012),  and  that  272 
these  interactions  are  driven  by  activity  of  local  popu-  273 
lations  of  neurons  in  the  cortex.  274 

275 

276 

4.  Discussion  277 

278 

This  study  shows  for  the  first  time  that  the  time-  279 
course  of  neural  activity  recorded  during  listening  to  2so 
music  is  conserved  across  individuals  in  cortical  ar-  2si 
eas  that  are  implicated  in  auditory  processing  —  not  282 
only  in  areas  that  process  low-level  aspects  of  audi-  283 
tory  stimuli  (i.e.,  areas  close  to  primary  auditory  cor-  284 
tex)  but  also  in  areas  that  process  higher-level  aspects  285 
of  auditory  stimuli  (i.e.,  auditory  associate  and  premo-  286 
tor  areas).  287 

The  results  shown  here  also  provide  exciting  new  288 
hints  about  fundamental  mechanisms  of  auditory  pro-  239 
cessing.  Previous  studies  have  provided  increasing  ev-  290 
idence  that  oscillations  in  the  alpha  band  modulate  291 
cortical  excitability  (Saalmann  et  al.,  2012).  Other  292 
studies  suggested  that  activity  in  the  high  gamma  band  293 
reflects  the  firing  rate  of  the  cortical  neuronal  popula-  294 
tion  beneath  each  electrode  (Miller  et  al.,  2009;  Miller,  295 
2010).  Based  on  these  existing  hypotheses  about  the 
mechanism  of  these  distinct  neurophysiological  pro¬ 
cesses,  one  may  expect  to  find  task-related  activity  in 
the  alpha  band  only  close  to  auditory  cortex,  and  task- 
related  activity  in  the  high  gamma  band  in  all  cortical 
areas  related  to  auditory  processing.  Our  results  pro¬ 
vide  the  first  direct  electrophysiological  confirmation 
of  this  expectation.  Specifically,  the  locations  of  sig¬ 


nificant  alpha-ISC  values  highlight  early  auditory  cor¬ 
tical  cortex  (Fig.  4c),  whereas  the  locations  of  signif¬ 
icant  gamma-ISC  values  highlight  cortical  areas  that 
have  been  proposed  to  be  activated  in  established  cor¬ 
tical  auditory  models  (Kaas  and  Hackett,  1999;  Za¬ 
torre  and  Belin,  2001;  Griffiths  and  Warren,  2002;  Za¬ 
torre  et  al.,  2002,  2007;  Hackett,  2008). 

Taken  together,  our  results  suggest  that  thalamo¬ 
cortical  projections  transmit  auditory  information 
from  the  thalamus  to  early  auditory  cortex,  from  where 
cortico-cortical  projections  relay  the  results  to  other 
(predominantly  peri- sylvian)  areas  that  process  more 
complex  auditory  features  (Zhang  et  al.,  2004;  Kumar 
et  al.,  2007;  Zatorre  et  al.,  2007)  (see  illustrative  yel¬ 
low  arrows  in  Fig.  4c).  Our  results  also  suggest  that 
stimulus-dependent  activity  in  early  auditory  areas  re¬ 
sults  in  delayed  feedback  to  the  thalamus.  In  this  view, 
the  thalamus  regulates  the  flow  of  auditory  informa¬ 
tion  into  the  cortex  based  on  this  feedback,  such  that 
increases  in  auditory  stimulus  intensity  result  in  sup¬ 
pression  of  alpha  activity,  which  in  turn  causes  an  en¬ 
hancement  of  cortical  excitability. 

These  results  have  several  important  implications. 
First,  our  approach  provides  a  novel  method  to  de¬ 
lineate  task-related  brain  activity,  which  is  currently 
evaluated  mostly  by  comparing  brain  activity  during 
a  task  to  brain  activity  during  rest  within  the  same 
subject.  Second,  the  approach  implemented  here  pro¬ 
vides  the  basis  for  functional  rather  than  anatomical 
co-registration  of  brains  of  different  individuals  of  the 
same  or  even  different  species  (Mantini  et  al.,  2012). 
Third,  the  finding  that  higher  stimulus  intensity  leads 
to  an  increase  in  cortical  excitability,  and  thus  to  a  pre¬ 
sumed  increase  in  the  consequences  of  the  stimulus, 
may  suggest  a  general  mechanism  to  prioritize  differ¬ 
ent  competing  sensory  streams,  and  could  thus  be  use¬ 
ful  in  explaining  different  types  of  stimulus -driven  or 
internally-generated  (covert  or  overt)  forms  of  atten¬ 
tion.  Finally,  the  concept  that  our  brains  speak  a  com¬ 
mon  “language”  when  they  are  coupled  to  the  external 
world  suggests  that  we  not  only  share  similarities  in 
our  bodies,  but  also  in  our  minds.  (Hasson,  2010). 
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Figure  1:  Subject- specific  brain  models  and  locations  of  implanted  electrodes.  The  brain  model  template  marked  with  a  star  shows  the 
approximate  location  of  implanted  electrodes  from  all  subjects. 
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Figure  2:  Intersubject  correlation  analysis.  (A)  Topographical  distribution  of  accumulated  negative  log  of  the  significance  of  the  ISC  values 
during  the  task  (left  panels)  and  rest  (right  panels)  conditions  for  ECoG  activity  in  the  alpha  (top  row)  and  high  gamma  (bottom  row)  bands. 
Values  larger  than  2  are  statistically  significant  at  a  confidence  level  of  99%  (see  vertical  line  in  color  bars).  (B)  Number  of  locations  with 
significant  ISC  values  at  each  frequency  band  and  for  task  (left)  and  rest  (right)  and  the  total  number  of  significant  locations  (277). 
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Figure  3:  Temporal  relationship  between  alpha,  high  gamma,  and  sound  intensity.  (A)  Representative  part  of  the  time  course  of  the  sound 
intensity  (black  trace)  and  the  average  time  course  of  ECoG  gamma  (blue  trace)  and  alpha  (red  trace)  activity.  (B)  Scatter  plots  illustrate  the 
relationship  between  these  three  variables  for  all  data  points,  as  well  as  regression  lines  and  their  95%  confidence  bands. 


Figure  4:  Relationship  between  anatomy  of  auditory  processing,  current  understanding  of  ECoG  physiology,  and  our  results.  (A)  Auditory 
pathways  between  the  medial  geniculate  nucleus  of  the  thalamus  and  primary  auditory  cortex.  Thalamic  neurons  (i.e.,  thalamo-cortical  (TC) 
and  reticular  (R)  neurons)  are  shown  in  gray,  auditory  cortex  (AC)  neurons  are  shown  in  blue,  and  their  interactions  are  shown  in  red.  (B) 
Top:  Exemplary  illustration  of  typical  power  spectral  density  of  ECoG  signals.  Activity  in  the  alpha  (8-12  Hz)  band,  indicated  by  the  red  band, 
reflects  thalamocortical  interactions.  Activity  in  the  gamma  (70-170  Hz)  band,  indicated  by  the  blue  band,  indexes  activity  of  local  populations 
of  neurons.  Bottom:  Average  time  course  of  the  first  120  sec  of  ECoG  activity.  The  blue  trace  gives  gamma  activity  —  the  red  trace  gives 
alpha  activity.  (C)  Outline  of  areas  with  significant  ISC  values  in  the  alpha  and  gamma  ECoG  amplitudes  given  by  the  results  in  Fig.  1.  Taken 
together,  our  results  are  consistent  with  the  hypothesis  that  auditory  information  reaches  the  cortex  in  areas  close  to  auditory  cortex,  from  where 
they  are  communicated  to  peri-sylvian  and  distinct  frontal  cortical  locations  (see  illustrative  yellow  arrows). 
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